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ABSTRACT 


The  purpose  of  this  document  is  to  provide  a  detailed  and  self-contained  deriva¬ 
tion  of  the  high  frequency  component  SLi  of  the  acoustic  surface  loss  formula 
for  a  rough  interface  of  impedance  mismatch.  We  thoroughly  examine  the  ex¬ 
isting  derivations  of  SLi,  point  at  their  flaws  and  suggest  the  corrected  form. 

Much  of  the  report  is  a  technical  review  in  which  we  critically  revisit  the 
fundamentals  behind  the  surface  loss  formula  and  cover  such  aspects  as  the 
basics  of  the  Kirchhoff  theory,  random  surfaces  and  wave  scattering  from  them, 
large  roughness  approximation  of  the  diffuse  field,  and  energy  flux  calculations. 
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Details  of  the  Derivation  of  the  Acoustic  Surface  Loss 

Formula 


EXECUTIVE  SUMMARY 

Tactical  decisions  in  ASW  often  have  to  be  based  on  projected  performance  of  sonars. 
Sonar  performance  prediction  in  turn  strongly  relies  on  acoustic  propagation  modelling. 
A  contemporary  model  for  underwater  sound  propagation  is  typically  a  complex  software 
which  implements  modules  describing  such  phenomena  as 

•  acoustic  propagation  in  an  ocean  with  a  general  sound  velocity  profile, 

•  reflection  of  sound  from  ocean’s  bottom  and  surface  (surface  and  bottom  loss  sub¬ 
models), 

•  volume,  bottom  and  surface  scattering  of  sound  and  the  resulting  reverberation. 

It  is  often  the  case  that  the  model  provides  more  than  one  sub-model  for  the  same  phe¬ 
nomenon  and  the  user  has  to  decide  which  of  them  is  most  appropriate  for  the  task.  Un¬ 
derstanding  the  fundamentals  upon  which  various  sub-models  are  based  helps  the  modeller 
to  make  more  informative  decisions  on  their  suitability. 

In  this  work  we  revisit  the  derivations  of  the  surface  loss  formula  used  by  the  “Beckmann- 
Spizzichino”  model  included  in  a  number  of  sonar  performance  prediction  tools.  Our  re¬ 
sultant  expression  for  the  high-frequency  surface  loss  term  SLi  is  in  variance  from  that 
currently  employed  by  the  model.  This  result  suggests  that  in  the  sonar  performance 
modelling  we  should  use  a  newly-obtained  expression  for  SLi  or  other  available  surface 
loss  models. 
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1  Introduction 


The  high  frequency  component  SLi  of  the  acoustic  surface  loss  formula  [2,  eqn.  3-8, 
p.  18]  broadly  used  in  acoustic  propagation  modelling  and  cited  in  the  literature  (e.g.  [1], 
[2],  [3],  [4]),  is  obtained  in  a  number  of  steps  and  has  a  deceptively  compact  and  simple 
form 


SLi  =  — 101og10(l  —  V), 


(1) 


where  V  is  given  by  (96)  and  depends  on  the  angle  of  the  incident  field,  the  surface 
correlation  length  and  the  RMS  value  of  the  surface  elevation  with  respect  to  the  mean 
plane. 

The  term  SLi  is  often  referred  to  as  the  “Beckmann-Spizzichino”  large-roughness  sur¬ 
face  loss,  presumably  because  some  elements  of  its  derivation  are  indeed  based  on  a  large¬ 
roughness  approximation  of  the  diffuse  component  of  the  reflected  field  obtained  in  [6, 
eqn.  47,  p.  87].  However,  the  authorship  of  SLi  can  hardly  be  attributed  to  the  authors 
of  [6].  Reports  [2,  1]  relate  the  origins  of  the  SLi-term  to  the  documentation  for  the 
RAYMODE  X  propagation  model,  but  provide  very  sketchy  demonstrations. 

Because  the  derivations  of  SLi  available  in  the  mainstream  literature  contain  errors, 
gaps  and  typographical  mistakes,  the  user’s  confidence  in  the  final  result  is  undermined. 
The  purpose  of  this  work  is  to  critically  revisit  the  fundamentals  behind  the  surface  loss 
formula  and  produce  a  document  with  its  detailed  and  self-contained  derivation. 

The  structure  of  this  document  is  as  follows: 

•  In  Section  2,  largely  based  on  [5],  we  review  the  Kirchhoff  theory  and  express  the 
complex  amplitude  of  the  scattered  field  in  terms  of  the  integral  over  the  mean  plane. 

•  In  Section  3  we  consider  the  coherent  and  diffuse  fields  and  the  associated  intensities 
and  examine  the  derivations  provided  in  [6]  and  [5]. 

•  In  Section  4  we  consider  the  asymptotic  behaviour  of  the  ensemble  average  of  the 
diffuse  intensity  at  large  values  of  the  roughness  parameter  and  confirm  the  key 
result  ([6,  eqn.  (47),  p.  87])  required  for  derivation  of  the  surface  loss  term  SLi.  We 
also  estimate  the  remainder  of  the  resulting  large-roughness  approximation  of  the 
diffuse  field. 

•  In  Section  5  we  carry  out  energy  flux  calculations.  We  estimate  the  total  flow  of 
the  diffuse  energy  and  use  the  result  to  derive  the  surface  loss  formula.  Step-by-step 
analysis  of  the  derivation  allows  us  to  identify  the  flaws  in  the  existing  version  of  the 
SLi  term  and  suggest  an  alternative  expression. 


2  Complex  amplitudes  of  scattered  fields 


The  complex  amplitude  of  an  acoustic  field  (scaled  pressure)  resulted  from  scattering 
of  an  incident  field  ^mc(r)  by  a  surface  So  enclosing  a  bounded  volume  satisfies  the  integral 
relation 
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where  if(r)  =  ^mc(r)  +  ipsc(r)  is  the  total  field,  ro  is  the  integration  variable,  dSo  is  the 
surface  area  element  at  ro,  d/dno  denotes  the  directional  derivative  along  the  normal  no 
to  So  at  ro  (see  Figure  1),  the  function  G( r)  is  given  by 


G(  r) 


eik\r\ 

47r|r|  ’ 


(3) 


and  k  is  the  wave  number.  Formula  (2)  is  referred  to  in  [5]  as  the  Helmholtz  integral  and 
its  derivation  is  based  on  the  Gauss  divergence  theorem  and  the  Green’s  formulae.  The 
following  incident  fields  are  typically  considered 


ipmc(r)  =  eiki-'r  (plane  wave  with  the  wave  vector  kjnc),  (4) 

?/>inc(r)  =  G( r  —  q)  (field  from  a  point  source  located  at  q),  (5) 

where  q  is  not  on  the  surface  Sq. 


Figure  1:  The  scattering  surface  So,  the  mean  plane,  the  radius-vector  tq,  and  vectors  kjnc 
and  no-  The  acoustic  field  is  considered  in  the  volume  zo  >  h(xo,yo),  where  zq  =  h(xo,yo) 
is  the  functional  form  of  So- 


The  proofs  of  the  equation  (2)  are  different  for  (4)  and  (5),  because  the  behaviour  of 
these  two  functions  is  different  -  the  field  (4)  is  regular  everywhere,  but  does  not  satisfy 
the  Sommerfeld  radiation  condition  [7,  eqn.  (9A-12),  p.  355],  whereas  the  field  (5)  satisfies 
the  Sommerfeld  radiation  condition,  but  is  singular  at  q.  We  revisit  the  derivation  of  the 
Helmholtz  integral  for  the  cases  (4)  and  (5)  in  Appendix  B. 

As  in  [5]  the  incident  field  is  assumed  to  be  the  plane  wave  (4).  The  Kirchhoff  approx¬ 
imation  adopts  the  assumption  that  the  values  and  normal  derivatives  of  the  incident  and 
reflected  fields  at  the  boundary  are  connected  in  the  same  way  as  in  the  case  of  reflection 
of  a  plane  wave  from  a  flat  surface 

^sc(r0)  =  7?o(ro)V’inc(ro)  and  =  -ifr( ro)^’J(ro)  ,  (6) 

(777-0  (777-0 

where  R( ro)  is  the  reflection  coefficient  at  ro-  Because,  for  the  plane  wave  (4),  we  also 
have 

difmc(r0)  ■ 

— - - =  i(kinc  •  nojV’  (r0)  , 

ono 
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we  obtain 

^(ro)  =  [l  +  JRo(ro)]V’inc(ro),  (7) 

=  i[l  -  i?o(ro)](kinc  •  n0)^inc(r0)  •  (8) 

on  o 

Application  of  (7)-(8)  implies  that  ro  is  not  shadowed  by  other  parts  of  the  scattering 
surface,  and  throughout  this  report  the  shadowing  effects  are  neglected.  Substitution  of 
(7)  and  (8)  into  (2)  expresses  ipsc(r)  in  terms  of  V,mc- 

Note  1.  The  Fresnel  reflection  coefficient  used  as  i?o(ro)  depends  on  the  angle  between 
k;nc  and  no  (see  Figure  1).  Since  no  depends  on  ro,  the  reflection  coefficient  also  depends 
on  ro,  even  if  the  reflecting  material  is  the  same  throughout  Sq.  However,  if  So  is  the  sea 
surface,  which  is  the  case  in  our  considerations,  the  Fresnel  reflection  coefficient  can  be 
assumed  to  be  —1. 


From  now  on  So  will  denote  a  patch  of  ensonified  surface.  This  patch  is  assumed 
to  be  centered  at  the  origin  and  its  size  d  =  maxroes0(?’o)  to  be  much  greater  than  the 
wavelength.  It  is  also  assumed  that  r>d  and  that  kd2 /r  <C  27t  (far-field  approximation). 
These  assumptions  allow  us  to  simplify  the  phase  in  G(r  —  ro) 

k\r  —  ro |  =  fc(r2  —  2r  •  rg  +  Vq)1^2  =  kr  —  ksc  •  ro  +  O(kro/r)  k,  kr  —  ksc  •  ro  ,  (9) 


where  ksc  =  kr/r,  and  to  reduce  G( r  —  ro)  and  its  normal  derivative  to 


G(r  -  r0)  «  - 


Jkr— iksc-ro 

47Tr 

gifc|r— r0| 


^G(r  -  ro)  =  -  4]f|r  _  ro| 

ifc|r— r0| 


ik  — 


1  \  <9|r  —  r0| 


47r|r  —  ro| 

gifcr— iksc-r0 

Attt 


ik  — 


r  —  r0|  y  dno 
1  \  (r  —  r0)  •  n0 


r  -  r0| 


r  -  r0 


-(n0  •  ksc)z . 


Substitute  (7),  (8),  (10)  and  (11)  into  (2)  to  obtain 

^SC(r)  =  '\..r  Js  [(^ok-  -  k+)  •  n0]  e*k  'r°dSo 


(10) 


(11) 


(12) 


where  k  =  kjnc  —  ksc  and  k+  =  k;nc  +  ksc. 


Let  xo,  yo  and  zo  be  Cartesian  coordinates  induced  by  the  mean  plane  as  shown  in 
Figure  1  and  let  ex,  ey  and  e2  be  the  unit  vectors  of  the  respective  axes. 


Figure  2  introduces  the  angles  6 1,  9-2  and  #3  which  describe  the  directions  associated 
with  the  incident  and  scattered  fields.  Note  that  the  vector  ksc  depends  on  the  radius- 
vector  r  of  the  observation  point. 

If  the  functional  form  of  the  surface  So  is  zo  =  h(xo,yo )  and  Sm  is  the  projection  of 
So  on  the  mean  plane,  then 


no 


— e 


dh 
x  dxo 


-  e  dh_  - l  e 
e Vdyo  +  “ 


(  dh 


2 


+ 


(13) 


UNCLASSIFIED 


3 


DSTO  -RR-0362 


UNCLASSIFIED 


zo 


Figure  2:  Angles  9\,  62  and  6*3. 


dSM  =  dS0(n0  •  e~) 


dS0 


the  integral  (12)  can  be  reduced  to  the  integral  over  Sm 

tp sc (r)  =  ^7  /  [(Rok~  -  k+)  •  ( -exhXo  -  eyhyQ  +  ez)]  elk  'r°dx0dy0 
•J  Sm 

Substitute 


kmc  =  A:(sin0iex  —  cos^ie^) 
ksc  =  k(sin  62  cos  9^ex  +  sin  62  sin  6^ey  +  cos  #2ez) 


into  (15)  to  obtain 


</>sc(r) 


ikeikr  [  (  dh  Bh  \ 

4vrr  JsM  \  dxo  +  dyo  7 

x  exp  {ik  [Ax 0  +  By0  +  Ch(x0,y0)}}  dx0dy0 


where 


A  =  sin  6\  —  sin  62  cos  6*3 
B  =  —  sin  6 2  sin  63 
C  =  —(cos  6\  +  cos  62) 

a  =  sin  6\ (1  —  Ro)  +  sin  62  cos  $3(1  +  Ro ) 
b  =  sin  02  sin  6*3(1  +  Ro) 
c  =  cos  $2(1  +  Ro)  —  cos  9\(1  —  Ro) 


(14) 


(15) 


(16) 

(17) 


(18) 


(19) 

(20) 
(21) 


(22) 

(23) 

(24) 


The  coefficients  A,  B  and  C  are  the  x,  y  and  z  components  of  the  vector  k  /k,  while  the 
coefficients  a,  b  and  c  are  the  x ,  y  and  z  components  of  the  vector  (k+  —  i?ok~)/A;. 
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Note  2.  The  derivations  above  are  no  more  than  a  conspectus  of  the  derivations  provided 
in  [5,  pp.  74-80]. 

Note  3.  Beckmann  &  Spizzichino  [6]  use  the  notation  v  =  k~  and  p  =  k+,  and  their 
formula  (18a)  on  p.  22  is  equivalent  to  (15). 

Note  4.  If  the  surface  elevation  h(x o,  yo)  =  0  for  all  xq  and  yo>  then  (18)  takes  the  form 


=  c  [  eikAx°+ikBy°dx0dy0  , 

4vrr  JsM 

and  if 

Sm  =  {(x0,yo)  |  -X  <x<X,  -Y  <y<Y}, 
and  Am  =  4 XY,  then 


jUpikr 

^oC(r)  =  V’otu  $1,  #2,  #3)  =  — T - c4m  sinc(fcAA)  sinc(£;Ay) , 

47T7’ 


(25) 

(26) 


(27) 


where  sinc(f)  =  sin(t)/f.  Equation  (27)  describes  scattering  of  a  plane  wave  from  a  flat 
2X-by-2Y  rectangular  plate. 

Note  5.  Beckmann  &  Spizzichino  [6]  use  the  normalised  amplitude  which  they  define  as 


?/>sc(r) 
$oc(r,9 1)  ' 


(28) 


where 

^  jUpikr 

%sc(r,0i)  =  ^cMi, 01,0)1^!  =  -^rAM  cos  0i.  (29) 

A  correction  factor  1/2  is  required  to  make  the  equation  (1)  on  p.  26  of  [6]  to  be  equivalent 
to  the  result  (18). 

Note  6.  It  is  noted  in  [5,  p.  76]  that  instead  of  equation  (2)  one  can  use  the  equation 


V,sc(r)  =  [ 

Js0  L 


<9?/>sc(ro) 
dun 


G(r-ro)-7T(ro) 


dG{ r  -  r0) 
dno 


dS0 


as  a  starting  point.  This  will  give 


(30) 


7/>sc(r)  =  /  i?0(k^  •  no)eik~-rodSo 

^  r  J  So 


instead  of  (12).  Equation  (18)  will  have  the  same  form,  but  for  the  coefficients  a,  b  and  c 
we  will  have  to  use  the  components  of  the  vector  (—Rok~/k). 


3  Coherent  and  diffuse  fields 

3.1  Coherent  and  diffuse  fields  and  their  intensities 

The  coherent  and  diffuse  fields  are  defined  as  follows 

Coherent  Field  =  (ipsc), 
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Diffuse  Field  =  ipsc  -  (V>sc>- 

Throughout,  (•}  is  the  mean  (expectation)  taken  over  all  permissible  realisations  of  the 
surface  Sq.  The  intensities  Ic  and  I d  of  the  coherent  and  diffuse  fields  are 

/C  =  |(V>SC>|2  and  Id=|V>sc-(^sc)|2. 

Using  these  definitions  we  obtain  the  relation  that  connects  (I)  [=the  mean  intensity  of 
the  total  field],  Ic  [=the  intensity  of  the  coherent  field],  and  (Id)  [=the  mean  intensity  of 
the  diffuse  field] 

(Id)  =  (l^sc|2)  -  l(^sc>|2  =  (-0  -  h  •  (31) 


Once  formula  (18)  has  been  obtained,  both  [5]  and  [6]  proceed  to  analysis  of  the 
composition  of  the  scattered  field.  They  consider  Sm  of  the  form  (26)  and  carry  out 
integration  by  parts  in  xo  and  yo  in  the  terms  containing  the  derivatives  dh/dxo  and 
dh/dyo  respectively,  as  a  result  of  which  ipsc(r)  can  be  represented  as 

^sc(r)  =  V’-e(r)  +  ^e(r) ,  (32) 

where 

=  -^2 F  [  eik^x^dx0dy0  (33) 

4nr  JsM 

+  ^  +  (34) 

<t>(x<h  yo)  =  Axq  +  By0  +  Ch(x0,  yo)  (35) 


and 


^e(r)  =  - 


ikeikr 
47T  r 


r  Y 


_  /  (  pikct>(X,y0)  _  pik<t>(-X,y0) 

kC 


ib_  rx 
+  kc  ]_x 


'-Y 


eik<j>(xo,Y )  _  ik^xo-Y) 


dyo 


dx  o 


(36) 


Note  7.  Equation  (34)  can  be  written  as 

k~  ■  (k+  —  i?0k_)  i?o|k^|2 


F  = 


2Ck2 


2Ck2  ’ 


(37) 


which  shows  that  the  term  'ip-e  is  invariant  with  respect  to  whether  equation  (2)  or  (30) 
is  used  as  a  starting  point  of  derivations  of  (34).  However,  due  to  different  definitions  of  a 
and  b  (see  Note  6),  the  term  if)e  does  depend  on  whether  (2)  or  (30)  are  chosen  for  deriving 
(36).  These  aspects  are  discussed  in  more  detail  in  [5,  p.  83]. 


Beckmann  &  Spizzichino  [6,  p.  31]  argue  that  the  term  ipe{ r)  can  be  neglected.  Ogilvy 
[5,  p.  86]  demonstrates  that  neglecting  the  term  '0e(r)  would  lead  to  an  erroneous  result 
when  evaluating  the  coherent  field.  Comments  in  [6,  p.  78]  near  equations  (25)  and  (26) 
confirm  that  Beckmann  &  Spizzichino  were  aware  of  the  necessity  to  account  for  the  “edge 
effects”  when  evaluating  the  coherent  field. 
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During  evaluation  of  the  mean  intensity  of  the  diffuse  field  both  [6]  and  [5]  assume 
that  the  edges  of  the  surface  h(xQ,yo )  are  fixed  at  the  mean  level,  that  is  h(±X,yo)  =  0 
and  h(xo,±Y)  =  0.  This  assumption  does  allow  us  to  neglect  the  terms  associated  with 
as  we  obtain 

(/d}  =  (|^-g|2)-|^-g)|2, 

however  it  invalidates  the  result  (-0SC)  =  ( e'lkCh)'i/jQC  (see  (4.34)  [5,  p.  86])  obtained  for  the 
coherent  field. 

In  order  to  evade  the  complications  associated  with  the  edge  term,  we  will  adopt  the 
Small  Slope  Approximation,  which  assumes  that  dh/dx o  and  dh/dyo  are  negligibly  small. 
Use  this  approximation  to  rewrite  (18)  as 


ibsc(r)  =  f  eik^Axo+Byo+Ch{xo'yo^dx0dy0 

4?rr  JsM 


(38) 


3.2  Assumptions  about  statistics  of  the  scattering  surface 


Suppose  that  the  random  surface  zo  =  h(xo,yo )  has  the  following  two-point  height 
probability  distribution 


p2(hi,h2,~R) 


2  <7g(l-Cg(J?))  / 


27raVl  -  Cl{R) 


(39) 


where  R  is  the  vector  drawn  from  the  first  to  the  second  point,  R  is  the  distance  between 
them,  h\  and  h2  are  random  elevations  at  the  first  and  the  second  point,  that  is 


R  —  Rx^x  Ry^yi  R  —  \J  4“  Ry , 

hi  =  h(x0,y0),  h2  =  h(x0  -  RXly0  -  Ry ). 


Following  [5,  6]  we  consider  Cq(R)  of  the  following  form 

C0(R)  =  e~R2/T2  .  (40) 


It  can  be  readily  verified  that  the  height  distribution  of  zq  =  h(xo,yo)  is  Gaussian 

P— h? /2a2 


/oo 

p2(h,  h2,H)dh2 

-OO 


a 


y/2^ 


(41) 


and  we  also  have:  (h)  =  0  and  (h2)  =  a2. 

It  can  also  be  shown  that  the  autocorrelation  coefficient  of  the  random  surface  zq  = 
h(xo,yo)  is  Cq(R),  that  is 


(h(x o,  yo)h(x0  -  Rx,  y0  -  Ry)) 
{h2(x0,y0)) 


C0(R)  ■ 


(42) 


Equations  (40)  and  (42)  allows  us  to  interpret  T  as  the  correlation  length. 
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The  characteristic  functions  of  p{h)  and  p2(h\ ,  h2,  R)  have  the  following  form 

/OO 

p{h)eishdh  =  exp  {-sV/2}  (43) 

-OO 

X2(si,s2lR)  =  (eisihl+is^)=  f  dh1dh2p2(huh2,  R)eis^+is^ 

J  R2 

=  exp  |~y(si  +  2C'0(jR)sis2  +  s|)|  (44) 

Formula  (44)  acknowledges  dependence  of  X'2  on  R  by  including  this  parameter  in  the  list 
of  its  arguments. 


3.3  Evaluation  of  intensities  and  their  means 


First,  applying  the  expectation  operator  to  (38),  we  obtain  the  coherent  field: 

<V>sc(r)>  =  -^c  f  eik^+By^x(kC)dx0dyo  =  ^c(r )X(kC) ,  (45) 

J  Sm 

where  V;oc(r)  defined  by  (25)  and  is  the  field  resulted  from  scattering  from  a  flat  surface 
{h  =  0). 

Note  8.  In  [5]  equation  (45)  is  obtained  without  resorting  to  the  small  slope  approximation 
(see  also  equation  (25),  [6,  p.  79]). 

The  intensity  of  the  coherent  field  thus  has  the  form 

Ic  =  \(^)\2  =  I0\X(kC)\2  =  I0e-9,  (46) 

where  g  is  the  roughness  defined  as 

g  =  k2C2a2  =  k2u2 {cos  6\  +  cos  02)2  ■  (47) 


Next  we  evaluate  the  mean  intensity  of  the  diffuse  field 
(Id)  =  (J)-/c 

1,2  2  p 

—  _  _  /  eik[A(xo-xi)+B(y0-y\)] 

(4yrr)2  JsMxSM 

X  ((jWho-h^  _  (jkCh^  (e-ikChoy^  dxodyodxidyi 

^.2C2  rX  j-Y  i-  2ir  /--Rmax(a"i,  yi,8) 

=  - ttt  /  dx  i  /  dyi  /  dd  dRR 

(4tt r)2  J_x  J-Y  Jo  Jo 

x  eJk [AR cos Q+br sin 9)  ^kC,  - kC ,  R)  -  x(kC)x{kC)\  (48) 

;.2C2  f'X  r-Y  rihnax  (an, 2/1,0)  _ 

=  -  2  /  dx i  /  dyi  /  dRRJo{kR\/ A2  +  B2) 

x  eJk [AR cos 6+br sin 0]  _  x(fcC')x(A:C')]  ,  (49) 


UNCLASSIFIED 


UNCLASSIFIED 


DSTO-RR-0362 


where  the  transition  from  (48)  to  (49)  is  based  on  the  formula 


1 

27T 


gi(x cos 0+y sin 0)^0  =  jq(^2^2)  . 


It  follows  from  (40)  that  the  function 

X2 (kC,  —kC,  R)  -  X(kC)x(kC)  =  e-k2c2F[i-c0(R)]  _  e~k2c2o2 


rapidly  decays  as  R  — ►  oo.  If  T  -C  VX2  +  Y1 .  then  we  can  simplify  (48)  by  replacing  the 
upper  limit  RmfLX(xi,  yi,9)  with  oo.  This  approximation  removes  dependence  on  x\  and  y\ 
from  the  integrand,  which  allows  us  to  carry  out  integration  in  those  variables  and  obtain 


(Id)  =  Am  I  "  Jo(kR\/ A2  +  B2) 


87rr2 


e~g[l-C0(R)]  _  e-g 


RdR 


(50) 


Note  9.  These  derivations  repeat,  with  minor  modifications,  the  derivations  provided  in 
[5,  p.  87]  and  [6,  p.  78-79].  Ogilvy  [5]  does  not  use  the  small  slope  approximation.  To  turn 
formula  (50)  into  the  result  obtained  in  [5]  we  simply  have  to  replace  c  with  2 F  (compare 
(33)  and  (38),  or  see  [5,  p.  82]),  which  gives 


k2F2 


(-^d) Ogilvy  =  27rr2  ^ M  /  Jo{kR\/ A2  +  B2) 


-g[l-C0{R)]  _  -g 


RdR 


(51) 


The  result  analogous  to  (51)  obtained  in  [6,  p.  79]  (equation  (31))  has  the  form 
2'KF2 

D{p}  =  J  MvxyT)[X2(vz,  - vz )  -  x(vz)x*(vz)\TdT  . 

This  has  to  be  corrected  first,  instead  of  A2  there  should  be  just  A.  Next  we  recall  that 
v  =  k_  =  (vx,vy,vz)  =  k(A,  B,C),  (see  [6]:  formula  (2)  p.  26),  so  vz  =  kC  and 

vxy  =  ^vl  +  v2  =  k\J  A2  +  B2.  (52) 

The  coefficient  T3  is  defined  by  equation  (3)  of  [6,  p.  29],  where  it  has  already  been  assumed 
that  the  reflection  coefficient  is  1.  Comparison  of  equation  (3)  of  [6,  p.  29]  with  equation 
(4.30)  of  [5,  p.  84]  shows  that 


F2s  =  F2  =  F2/cos261, 


where  Fq  =  F  and  F  is  given  by  (34)  (or  (37))  in  which  Rq  =  —1,  that  is 


Fq  =  F  = 


1  +  cos  6 1  cos  02  —  sin  0\  sin  62  cos  63 
cos  9\  +  cos  02 


(53) 


(54) 


If  we  divide  (51)  by  fyosc|2  (see  (29))  and  take  the  above  considerations  into  account,  we 
readily  obtain  that  (51)  and  the  result  for  D{p}  obtained  in  [6,  p.  79]  and  reproduced 
above  are  equivalent. 
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Equations  (50)  and  (51)  become  the  same  in  the  specular  direction,  for  which  we  have 
A  =  B  =  0,  F  =  c/2  =  Rq  cos  0±  =  —  cos  9\,  g  =  go  =  Ak2a 2  cos2  9\  and  so 


/T  \  it  \  k2  cos2  6*i  _on  f 

(Id)  =  (^d)spec  =  - 2tTT^ - ^A/e  9  J 


Using  the  expansion 


exp  g0  e 


0-R2/t2\ 


egoCo(R)  _  ^ 

00  n  -nR2/T2 


RdR . 


E 

71=0 


%e 


n\ 


equation  (55)  can  be  transformed  to 

k 2  cos2  9i 

(f  d)  spec  —  2/7TT2  AmC 


-go 


fOO  oo  n c—nR?  IT2 

'  ^  yoe 


E 


RdR 


T2k 2— 2 


COS^  (71 


71=1 
OO 


ni 


47T7-2 


4/e-®  Y, 


71=1 


go 

n!n 


which  gives  the  following  expression  for  the  mean  square  normalised  amplitude 
„  1 2 \  |4sc(r,0i)| 2e-so  +  (/d)sre,n  T2tt 


/spec 


l^osc(u  ^Ol2 


/sPec  =  g— 90  + 


a— 90 


E 

71=1 


go 

n!n 


(55) 


(56) 


(57) 


Medwin  [8]  studied  scattering  in  the  specular  direction  and  equation  (57)  has  the  same  form 
as  his  key  equation  (8)  (see  [8,  p.  8]).  Note  that  Medwin’s  derivation  of  (57)  is  incomplete 
and  contains  errors1.  Equation  (57)  can  also  be  readily  obtained  from  equation  (57)  of 
[6,  p.  88]  if  we  observe  that  in  the  specular  direction  g  =  go,  vxy  =  0,  Eg s  =  1  and 
po  =  sine ( Ur X)  sinc(u2/y)  =  1. 

Note  10.  The  statements  (4.56)  of  [5,  p.  92]  and  (44)  of  [6,  p.  87]  are  erroneous.  They 
both  interpret 

Jo(kRy/ A2  +  B2)e~9[1-Co{R)]RdR  (58) 

as  (I)  when  g  — >  oo.  The  integral  in  (58)  is  divergent,  so  (58)  is  meaningless  and  can  be 
interpreted  as  neither  (I),  nor  (Id)-  Similarly,  the  expression 

k2F2 
2tt  r2 

is  equally  meaningless  and  cannot  be  interpreted  as  Ic.  Also  note  that,  because  the 
integrals  in  (58)  and  (59)  are  divergent,  we  cannot  replace  in  (48)  the  upper  integration 
limit  Rmax(xi,yi,9)  with  oo  if  (I)  and  Ic  are  considered  separately. 


/•uu 

Am  /  J0(kR\J A2  +  B2)e~9 RdR 


(59) 


k2F 2  r00 

^Am  L 


4  Asymptotic  behaviour  of  (Id) 

In  this  section  we  verify  the  asymptotic  formula  describing  the  behaviour  of  (Id)  when 
the  roughness  parameters  g  =  k2C2a 2  is  large  for  “gently  undulated”  surfaces  [6,  p.  194], 
which  satisfy  Aa2 /T2  <C  1. 

1The  last  equation  on  page  7  of  [8]  before  formula  (8)  is  stated  on  page  8  is  incorrect. 
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The  main  idea  of  the  proof  is  outlined  in  [6,  p.  87]  and  consists  in  replacing  the  integral 
in  (50)  with 

J0(kR^A2  +  B2)e~9R2/T2RdR  .  (60) 

Beckmann  &  Spizzichino  argue  that  the  main  contribution  to  the  integral  in  (50)  is  due  to 
integration  in  the  neighbourhood  of  R  =  0,  where  Cq(R)  can  be  approximated  by  ~  1  —  ^ ■ 
Using  this  approximation  and  dropping  and  adding  terms  insignificant  when  g  — >  oo  leads 
to  the  expression  (60),  which  can  be  evaluated  analytically.  The  resulting  estimate  for  the 
normalised  intensity  is  equation  (47)  of  [6,  p.  87]  reproduced  below 


(pp*) 


7 tF'2T2 

Ag 


exp 


vlyT2\ 

4  g  ) 


(eqn.  (47),  [6,  p.  87]) 


(*) 


where  we  have  to  recall  that  vxy  is  given  by  (52)  and  F  must  be  interpreted  as  Tbs  (see 
(53)).  Also,  A  in  (*)  must  be  interpreted  as  Am-,  that  is  the  area  of  Sm- 

Note  11.  If  we  had  k2T2/4g  <C  1,  then  retaining  the  asymptotic  formula  in  the  form 
(*)  would  be  meaningless,  as  the  principal  term  would  be  simply  the  coefficient  before  the 
exponential.  In  our  case  k2T2/4g  =  T2 /(Aa2C2)  isn’t  a  small  value,  so  (*)  is  exponentially 
small  unless  A2+B2  is  equal  to  or  close  to  zero,  which  corresponds  to  the  specular  direction. 


Consider  the  integral  in  (50) 


S(g)  = 


e-9(l-exp(-T2/T2))  _  e-g 


rdr 


(61) 


where  we  have  switched  from  R  and  k\J A2  +  B2  to  more  concise  r  and  vxy. 

To  verify  the  estimate  (*)  and  obtain  the  asymptotic  formula  for  (Id)  we  have  to  study 
the  behaviour  of  the  integral  (61)  at  large  g.  This  is  done  in  Appendix  C  where  the 
following  result  is  verified  (see  equation  (C14)) 


T2 

S(g)  =  7T 
2  g 


exp 


as  g  oo. 

Equation  (50)  and  estimate  (62)  give 
,T  .  k2c2T2 

(Id)  =  V7- - 


16irr2g 


vlyT 2 
4  g 


+  0(g 


-i\ 


exp 


VZXyT 


4  g 


+  0(g 


(62) 


(63) 


Note  12.  We  have  been  using  the  small  slope  approximation.  If  this  approximation  had 
not  been  used,  the  principal  term  of  the  asymptotic  formula  for  (Id)  would  have  the  form 
(see  [5,  p.  92],  equations  (4.57)  and  (4.58)): 


/rx  k2F2T2  A  (  v2xyT2\ 

(4d)ogilvy  —  47rr2g  Am  eXP  y  4 g  J 


(64) 


Equation  (64)  can  be  obtained  from  the  principal  term  of  (63)  by  simply  replacing  c  with 
2 F  (compare  (33)  and  (38),  or  see  [5,  p.  82]). 
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Note  13.  If  we  divide  (64)  by 

~  k2  A2 

l^oSC|2  =  77, — ¥2  c°s2  Oi  (see  (29)) 

( 2-Kry 

and  use  (53)  we  obtain  equation  (47)  of  [6,  p.  87]  (see  (*)  on  page  11). 

Note  14.  The  mean  power  reflection  coefficient  is  defined  as  (|'i/,sc|2)/'ko>  where  S&o  is  the 
amplitude  of  the  incident  wave  [6,  p.  89].  In  our  case  \ko  =  1)  so  (I)  (~  (Id)  when  g  1) 
can  be  also  interpreted  as  the  mean  power  reflection  coefficient.  Beckmann  &  Spizzichino 
(see  [6,  p.  89],  equation  (63))  provide  the  following  result  for  the  principal  term  of  the 
mean  power  reflection  coefficient: 

<|^i|2>  =  ^-Tcot2  A)exp  ^-tan2/goJ  (ecln-  (63),  [6,  p.  89])  (f) 


where  the  angles  (3q  and  (3  are  such  that 


2(7  Vtii 

tan  /3o  =  —  and  tan  (3  =  — -  . 
T  vz 


(65) 


If  we  correct  (f )  by  dividing  its  right-hand  side  by  4  and  use  (65)  together  with  the  relations 


F  =  — 


Po|k 


— 12 


Rqv 


(as  Rq  =  —1) 


2  Ck2  2  vzk  2  vzk 

(see  (37))  and  g  =  a2C2k2  =  a2v2,  we  obtain  that  (f)  and  (64)  are  equivalent. 
Note  15.  Because  we  have 


2,2  2 

^ xy  4“  Vz  V 


=F  —  —  — -M-  +  1  —  tanz  (3+1  — 


1 


cos2  (3 


we  can  use  an  alternative  form  of  (64) 


Am  cot2  (3o 

(  d)°gilvy  “  4vrr.2  cos4  p  eXP 


tan2  (3 
tan2  (3q 


(66) 


This  would  coincide  with  equation  (7)  of  [6,  p.  251]  if  we  take  into  account  that  there  is  a 
missing  factor  1/4. 

Note  16.  While  (3  can  be  readily  interpreted  as  the  angle  between  the  z-axis  and  the  vector 
k“  =  kinc  — ksc,  interpretation  of  (3q  is  less  straightforward.  Parameter  tan2  (3q  is  connected 
with  the  mean  square  slope,  which  is  defined  for  a  two-dimensional  surface  zq  =  h(xo,yo) 
as  the  ensemble  average  of  the  square  tangent  of  the  angle  between  the  normal  no  to  the 
surface  So  and  the  vector  e2.  It  can  be  readily  shown  that  the  mean  square  slope  can  be 
written  as  (h20  +  h20),  where,  as  before,  the  subscripts  denote  the  partial  derivatives  in 
respective  variables  and  (•)  denotes  the  expectation  operator. 

To  establish  the  connection  between  ( h20+h20 )  and  tan  (3q  consider  the  spatial  processes 
hxo(xo^yo)  and  hyo(xo,yo)  derived  from  the  random  process  h(xo,yo). 

Let  Po(k),  P\  (k)  and  P2(k)  be  the  power  spectrums  of  the  processes  h,  hXQ  and  hyo, 
and  let  iio(R),  A'i(R)  and  /^(R)  be  their  autocorrelations.  By  (40)  and  (42),  ATo(R)  is 
a2  exp(— R2/T2).  From  the  Wiener-Khintchine  relations,  Po(k)  has  the  following  form 

Po(k)  =  —^2  r  dRx  r  dRyKoiny™  =  4 

)  J  —  oo  oo 
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Using  the  properties  of  the  Fourier  transform  and  the  definition  of  the  power  spectrum, 
we  obtain  -Pi(k)  =  fc2.Po(k)  and  P2(k)  =  fc2Po(k).  By  the  Wiener-Khintchine  relations  we 
have 

/oo  /*oo 

dkx  /  dkyP£(k)eikn  (£  =1,2). 

-OO  J  —  OO 


Hence 


and  so 


K1(R)  =  --^Ko(R)  and  K2(R)  =  ~-^K0(R) 

r)2  2rr2 

(hlQ)  =  Ad(0)  =  -  ^o(R)  =  w 

R— 0 

82  9  rr2 

(^o)  s  x2(o)  =  -  w^o(R)  =  • 

y  r=o 


This  gives  us  the  equation  connecting  the  mean  square  slope  of  the 
tan2  fa 


{h2Xo  +  h 


4  cr2  9 
=  -f2  =  tan-  po- 


scattering  surface  and 


Using  optical  measurements  Cox  &  Munk  [9]  obtain  regression  formulae  expressing  (/i20), 
K)  and  (hlo  +  hlo)  in  terms  of  the  wind  speed.  This  was  used  by  Medwin  [8]  to  generate 
correlation  lengths  T  for  the  wind  speeds  of  interest.  Note  that  Medwin  [8]  appears  to 
erroneously  apply  the  Cox  &  Munk  [9]  regression  formula  for  {h^.Q+hy0)  to  {hXo}  =  2a2 IT2. 


5  Energy  flux  calculations 

In  this  section  we  examine  the  rates  of  outgoing  energies  of  the  diffuse  and  coher¬ 
ent  components  of  the  scattered  field.  First,  in  Section  5.1,  we  provide  calculations  for 
scattering  by  a  flat  horizontal  plate.  These  calculations  show  that  all  energy  incident  on 
a  perfectly  reflecting  flat  plate  is  entirely  reflected  in  a  narrow  beam  directed  along  the 
specular  axis.  Next,  in  Section  5.2,  we  consider  a  patch  of  a  general  gently  undulated 
surface  and  evaluate  the  total  flux  ntj  of  the  diffuse  energy  of  the  scattered  field  for  the 
situation  when  the  roughness  parameter  is  large.  We  show  that  if  the  incident  wave  vector 
k;nc  is  separated  from  the  horizontal  direction  then  the  total  diffuse  flow  is  equal  to 
the  incident  energy  flow  no- 

For  wave  vectors  k;nc  closely  aligned  with  the  horizontal  direction,  the  expression  for 
Ha  is  more  complicated  but  provides  us  with  an  explanation  as  to  why  the  term  V  in 
the  high-frequency  surface  loss  formula  has  the  form  (96).  The  procedure  described  in  [2, 
p.  54-56]  uses  the  total  diffuse  flow  in  order  to  obtain  the  coherent  reflection  coefficient.  If 
we  strictly  follow  this  procedure  we  obtain  the  required  expression  for  SLi.  However,  we 
believe  that  equation  (3B-2)  of  [2,  p.  54],  which,  in  the  notation  of  [2],  has  the  following 
form 

/*27T  /*tt/2 

R+  /  <7 cos  (j}rd(f>rd9r  =  1,  (3B-2,  [2,  p.  54]) 

Jo  Jo 
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is  valid  only  for  a  vertically  incident  wave.  For  arbitrary  grazing  angles  fa  equation  (3B-2, 
[2,  p.  54])  should  read 

re2'K  n'K  1 2 

R  sin</>i+  /  /  <7  cos  fadfadd?  =  sin  fa  . 

Jo  Jo 

The  new  expressions  for  SLi  resulted  from  this  correction  are  provided  in  Section  5.3  in 
Notes  18  and  19. 


5.1  Energy  flux  calculation  for  reflection  from  a  flat  plate 

The  field  due  to  reflection  from  a  rectangular  flat  horizontal  plate  (26)  is  given  by  (26) 
(see  Note  3) 

jbpikr 

'0nc(r)  = - - - cAm  sinc(AuLA)  sinc(/cAT)  . 

47t  r 

The  total  rate  of  energy  reflected  from  the  plate  (26)  is  given  by 

pn  /-7r/2 


rin  = 


d03  J  d02  sin  92r2\'iplc\2 

de3 £'2 de2Sme2?Ak  ( ( sin2(fct?) 


nkXA 2  )  V  7 rkYB2  ) 


(67) 


To  evaluate  the  integral  (67)  note  that  both  kX  and  kY  are  large  and  that  the  sequence 

sin2(?nf) 


Sm{x)  = 


nmt 2 


(68) 


is  a  (^-sequence  as  m  — >  oo,  that  is 


lim 

m— kx) 


Sm(t)f(t)dt  =  /( 0) 


(69) 


for  any  smooth  f(t)  vanishing  outside  a  finite  segment.  The  result  (69)  can  be  readily 
verified  if  we  make  use  of  the  following  facts: 


lim  Sm(t)  = 


oo,  if  t  =  0  (as  5m(0)  =  m/n) 
0,  if  t  ±  0 


and 


sin2  mt 
irmt2 


dt  =  1 . 


(70) 


(71) 


By  (70),  the  integrand  of  (67)  vanishes  outside  a  small  neighbourhood  of  the  point  where 
A2  +  B2  =  0  (that  is,  where  we  have  9o  =  9\  and  93  =  0),  which  corresponds  to  the 
specular  direction2.  Therefore  the  integration  domain  can  be  reduced  to  a  small  cone  of 
directions  which  surrounds  the  specular  axis.  In  order  to  make  use  of  (69)  we  change  from 
02  and  #3  to  new  integration  variables 


d  =  A(02,93),  v  =  B(9  2,9s): 


(72) 


"This  fact  and  equation  (46)  explain  why  the  coherent  field  is  described  as  a  specular  spike  [10] 
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where  A  and  B  are  given  by  (19)  and  (20). 

The  Jacobian  of  the  transformation  defined  by  (72)  is 


J(0  2,  Os) 


9{£,v) 

d(02,03) 


"402 

"4e»3 

—  cos  02  cos  03  sin  0-2  sin  03 

Be2 

Bb3 

—  cos  02  sin  03  —  sin  02  cos  03 

and  we  have 


det  J  =  cos  02  sin  62 


1  _  1 
detJ|  cos  02  sin  02 


Equation  (67)  transforms  to 


nn 


dij 


c2  Am 

4  COS  02  62=62^,1) 
63=63(1,1) 


( sin 2(kX£)  \ 

v  vr  kxe  ) 


(  sin  2(kYrj)\ 
\  irkYrj2  J 


(73) 


where  integration  is  carried  out  in  small  neighbourhoods  of  £  =  0  and  7]  =  0.  The  pair 
02  =  02(C)7?)  and  03  =  03(C)7?)  used  in  (73)  is  a  transformation  inverse  to  (72). 

Now  use  (70)  to  obtain 


n 


0  — 


c2Am 
4  COS  02 


62=62  (0,0) 
03=03(0,0) 


Am  cos  0i , 


(74) 


where  we  have  used  (24)  (with  the  substitution  Rq  =  —1)  and  the  fact  that  02(0,0)  =  0i 
and  03(0,  0)  =  0. 

The  incident  field  elkinc'r  results  in  acoustic  energy  flowing  in  the  direction  k;nc.  The 
rate  at  which  this  flow  delivers  the  energy  to  the  horizontal  plate  of  area  Am  is  Am  cos  0i. 
The  result  (74)  is  intuitively  satisfying  as  it  shows  that  under  the  Kirchhoff  approximation, 
upon  which  the  solution  ipff  is  based,  all  energy  delivered  to  a  flat  perfectly  reflecting  plate 
is  entirely  reflected  in  a  narrow  beam  of  directions  centred  at  the  specular  axis. 


5.2  Calculation  of  the  total  diffuse  energy  flux 


The  rate  at  which  the  surface  So  scatters  the  incident  field  e*kinc’r  into  the  diffuse 
energy  is  given  by  the  equation 


nd  =  /  d03  f  d02  sin  02r2(Id) 

J-ir  Jo 


(75) 


where  (Id)  is  given  by  (50).  If  the  roughness  parameter  g  is  large,  equation  (63)  (or  (64), 
if  we  do  not  use  the  small  slope  approximation)  is  valid  and  (75)  takes  the  form 


nd  =  Am  cot-  p0  / 

J  —  7T 

where  Q{02,0s)  is  given  by 

k2c2T2  sin  02  4ct2 


r/z  f  tan2  B 

d03l  ‘is 


Q(02)  0 3)  — 


4  cos2  0i  sin  02 


cos2  0i  sin  02 


167T<7  T 2  47r(cOS  01  +  COS  02)2  7r(cOS  01  +  COS  02)2 


(76) 


(77) 
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in  the  small  slope  approximation,  and 


Q(6  2, 63) 


k2F2T 2  sin  02  4cr2  _  F2  sin  02 

Ang  T 2  7r(cos  9\  +  cos  6*2) 2 

(1  +  cos  0\  cos  d2  —  sin  0\  sin  02  cos  6*3) 2  sin  02 
7r(cos  6\  +  cos  02) 4 


if  we  use  (64)  obtained  without  assuming  the  small  slope  approximation. 
The  integral  (76)  can  be  estimated  at  large  values  of 

T 2 

g  =  cot2  /3o  =  ^2 

using  a  multi-dimensional  form  of  the  Laplace  method  [11,  Ch.  2].  Let 

$(02,03)  =  —  tan2  j3  . 


(78) 


(79) 


(80) 


The  maximum  of  value  of  $(02,03)  on  11  =  {(02,03)  £  [— vr,vr]  x  [0, 7r/2] }  is  zero  and  it  is 
achieved  if  and  only  if  both  A  and  B  are  zeros  (see  (65)  and  (19)— (21)),  which  happens  at 
the  point  (62,83)  =  (0i,O)  corresponding  to  the  specular  direction. 

If  6\  is  separated  from  7t/2,  then  (0i,O)  is  an  internal  point  of  the  integration  domain 
and  application  of  the  Laplace  method  (e.g.  [11,  Theorem  4.1,  p.  74])  gives 

■^,^(*,0)}  (7)  V|  det  '  <81) 


where  g  and  are  given  by  (79)  and  (80),  and  H(62, 63)  is  the  Hessian  matrix  of  $(#3, 63) 
defined  as 

$61202  $0203 

$0203  $0303. 

It  can  be  shown  that 


H  (62,63)  = 


H(61,0) 


0 

—  \  tan2  6\ 


det  i7(0i,O) 


1  2 

-  tan 
4 


6 1. 


(82) 


Substitution  of  d2  =  0i  and  83  =  0  into  (77)  and  (78)  gives 


Q(8 1,0) 


sin  0i 
4-7T 


(83) 


for  both  cases.  Using  (82),  (83)  and  the  fact  that  <h(0i,O)  =  0  we  obtain  that  (81)  has  the 
form 


27tQ(0i,O)Am 
VI  det  i7(0i,  0)| 


Am  cos  0i  . 


(84) 


If  0i  is  allowed  to  approach  n/2  we  can  use  the  following  argument.  Since  in  (76) 
tan-2  /?o  is  large,  the  main  contribution  to  the  integral  is  made  by  permissible  directions 
which  belong  to  a  small  cone  K.5  around  the  specular  axis  (02,  03)  =  (0i,  0)  where  tan2  [3  =  0 


JCs  —  {(62, 83)  £  H  |  0i  —  6  <  02  <  0i  +  6,  —5  <  63  <  8} 
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Because  (6*1,0)  is  a  stationary  point  of  tan2/?  and  (cf.  (82)) 


i 52  ttan2  P\ 


02=01  2 
03=0 


— 2  [tan2  /?]  =-tan20! 


02=01  2 
03=o 


302  6*3 


[tan2  /? 


we  can  use  Taylor’s  expansion  and  replace  tan2  /?  inside  K,$  with 

tan2  ^  [(02  —  0i)2  +  tan2  0i0|]  (88) 

Approximating  <5(02,03)  by  its  value  at  (02,03)  =  (0i,O)  (see  (83))  and  using  (88)  we 
transform  (76)  to  the  form 

r*  r . 

4tt  y_,5 

The  presence  of  the  large  parameter  in  the  rapidly  decaying  exponential  allows  us  to  extend 
the  limits  of  integration  without  introducing  a  considerable  error 


nd  ^ 


Am  sin  0i  2 


cot2/?0  [°°  d03  ^  d02e-C-^{(e2-9l)2+tan2ei93'\ 

0-00  J—  oo 


Changing  in  (89)  to  new  variables 


£  =  7,  cot(/?o)(02  -  0i),  0  =  ^  cot(Aj)  tan(0i)03 


we  obtain 


nd  ~ 


Am  cos  0i 


Am  cos  0i 


e  T*  d?? 


•  i  (tt/2 — 0i )  cot /3q 


e  ^  d£ 


e  ^  d( 


J - j(7t/2—  0l)  cot  00 

=  ^  Am  cos  01  erfc  ^  -  0i)  cot  /30^ 


where  erfc(rc)  is  the  complementary  error  function 


2  r°°  2 

erfc(tc)  =  —j=  /  eTl  dt . 
V7*-  J w 


Equation  (90)  takes  a  more  compact  form  if  it  is  written  in  terms  the  grazing  angle 
a\  =  7t/2  —  0i 


nd  ~  -Am  sin  a \  erfc  — 


aq  cot  fa 
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Note  17.  Let  us  summarise  the  assumptions  used  in  the  derivation  of  (91).  In  addition  to 
assumptions  that  guarantee  applicability  of  the  Kirchhoff  theory  (see  [5,  Section  4.2])  we 
have  also  assumed  the  following: 

1.  The  two-point  height  distribution  of  the  scattering  surface  is  Gaussian  (see  (39))  and 
the  autocorrelation  coefficient  has  an  exponential  form  (40). 

2.  The  correlation  distance  T  is  much  smaller  than  the  size  of  the  patch:  T  <S  X 
and  T<  7.  This  is  required  for  stochastic  treatment  utilising  the  stationarity  of 
the  random  surface  [5,  pp.  24-26]  and  this  assumption  was  also  used  when  (49)  was 
transformed  to  (50). 

3.  The  value  of  the  roughness  parameter  is  large:  g  =  Ak2cr2  cos2  9\  1  (here  62  =  0\). 

4.  We  consider  two-dimensional  gently  undulated  surfaces,  for  which  tan2  L 

Note  that  by  Note  16  this  requirement  is  consistent  with  the  adopted  small  slope 
approximation. 

5.  Shadowing  effects  are  neglected. 


Numerically  the  complementary  error  function  is  efficiently  evaluated  using  the  Cheby- 
shev  fitting  (e.g.  [12,  pp.  220-221]),  so  the  form  (91)  hardly  required  any  further  reduction. 
However,  the  following  asymptotic  formula3  for  erfc 

e-n2  / 1  \ 

erfc(— u)  =  2 - —  (  — b  0(u~ 2)  )  ,  as  u  — »  00,  (92) 

W  ) 


has  resulted  in  further  transformation  of  (91)  to 


nd  ~  Am  sin  aq 


1  - 


exp 


cot2  (3q  a2 
4 


7 r  aq  cot  fa 


Formula  (92)  is  valid  when  u  is  large;  when  u  — >  0,  the  right-hand  side  of  (92)  tends  to 
—00,  whereas  we  must  have 

erfc(— u)  >  erfc(0)  =  1 .  (93) 

The  fact  that  (92)  is  singular  at  u  =  0  and  the  requirement  (93)  have  apparently  resulted 
in  the  following  formula  (e.g.  [1,  2]) 


n d  —  Am  sin  a\  max 


(  erfc(0) 


1  - 


exp 


cot2  /?0  cuf  ^  n 
4 


7T  oq  cot  /3q 


.  sin  aq 

~  Am  max  |  — - — ,  smaq 


1  - 


exp 


cot2  (3q  af 
4 


7T  aq  cot  /?o 


(94) 


^Formula  (92)  can  be  readily  obtained  using  integration  by  parts  (e.g.  [13,  Ch.  3,  Section  1]). 
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5.3  The  acoustic  surface  loss  formula 


The  high  frequency  acoustic  surface  loss  formula  used  in  propagation  modelling  (e.g. 
[1,  p.  45],  [2,  p.  56],  [3])  has  the  following  form 


where 


SL,  =  — 101og10(l  —  V) 


V  =  max 


sin  aq 


sin  aq 


1  - 


exp 


cot2  /3o  ctf  ^  n 
4 


7 T  a  1  cot  Po 


(95) 


(96) 


The  derivations  of  equation  (94)  clarify  the  nature  and  the  composition  (equation  (96))  of 
the  term  V  that  appears  in  (95). 

Note  18.  Evaluation  of  the  surface  loss  SLq  using  the  definition  (cf.  [14,  eqn.  (9.6.1), 
p.  199]  or  [15,  eqn.  (4),  p.  230]) 


SLq  =  — 101og10  Hlh—Ell ,  (97) 

Ho 

where  IIo  and  Ilq  are  given  by  (74)  and  (94),  gives  a  result  different  from  (95).  Namely, 
we  obtain 

SLi  =  —  101og10(l  —  U/ sin  aq)  (98) 

with  V  having  the  same  form  (96)  as  in  (95). 

Note  19.  If  in  (97)  we  use  for  Ilq  equation  (91)  instead  of  (94),  we  obtain 


SLi  =  —  101og10 


-  erfc 
2 


oq  cot  (3q 


(99) 


6  Concluding  Remarks 

The  findings  of  Section  5  demonstrate  that  the  currently  used  form  of  the  surface  loss 
term  SLi  is  erroneous.  Our  analysis  identifies  at  least  two  flaws  in  the  derivations  leading 
to  equations  (95)  and  (96).  One  of  these  flaws  is  an  apparent  fix  (94)  required  as  a  result 
of  the  application  of  asymptotic  expansion  (92)  outside  its  validity  region.  The  second 
flaw  is  more  fundamental  and  consists  of  an  incorrect  assumption  that  the  energy  flux 
balance  equation  for  a  wave  of  vertical  incidence  can  be  applied  to  an  arbitrarily  incident 
wave.  Our  corrected  expression  for  SLi  is  given  by  equation  (99). 

Equation  (99)  follows  from  equation  (97)  which  implies  that  the  coherent  energy  flow 
can  be  evaluated  from  the  total  and  diffuse  energy  flows.  If  we  use  an  alternative  definition 

SLC  =  —10  log10 

Ho 


where 


nc 


f  71-/2 


dd 2  sin  02  rlc 


(100) 
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Frequency  4000  Hz 


Grazing  angle  a  (deg) 


Frequency  8000  Hz 


Figure  3:  Comparison  of  SLC  (eqn.  (101)),  SL  (eqn.  (102)),  and  the  different  forms 
of  SLi  given  by  (95)-(96),  (98),  (99)  and  (101)  for  two  wave  frequencies  4  kHz  (a)  and 
8  kHz  (b). 
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and  Ic  is  the  intensity  of  the  coherent  field,  then  equation  (46)  and  an  argument  similar 
to  that  used  in  Section  5.1  give 


SLC  =  10  log10  e9  .  (101) 

This  result  is  attributed  to  Eckart  (see  [16,  eqn.  (26),  p.  569]  and  [3]). 

Formula  (63)  shows  that  for  gently  undulated  surfaces  the  diffuse  field  rapidly  vanishes 
as  we  move  from  the  specular  direction.  This  is  somewhat  in  contrast  with  the  view  (e.g. 
[5,  p.  5-6])  that  the  specular  field  is  coherent  and  the  diffuse  field  has  a  wide  angular  spread. 
Medwin  [8]  used  [6,  eqn.  (57),  p.  88]  to  show  that  at  small  values  of  roughness  parameter 
g  the  coherent  component  dominates  in  the  specular  direction,  whereas  at  large  values  of 
g  the  field  in  the  specular  direction  is  dominated  by  the  diffuse  component.  If  the  acoustic 
modelling  uses  incoherent  addition  of  eigenray  contributions,  then  the  reflection  coefficient 
considered  in  [8]  (see  equation  [8,  eqn.  (8)]  or  our  equation  (57))  and  the  associated  surface 
loss 

SL  =  —10 log10  (J)l^=e^e3=o  =  _ioioglo(|p|2)spec  (102) 

h 

should  be  used  instead  of  the  coherent  surface  loss  (101)  or  its  alternatives,  such  as  (98) 
or  (99).  Note  that  the  area  of  the  scattering  patch  is  among  the  inputs  of  (57). 

Figure  3  compares  SLC  (eqn.  (101)),  SL  (eqn.  (102)),  and  the  different  forms  of  SLi 
given  by  (95)— (96) ,  (98),  (99)  and  (101)  for  two  wave  frequencies  4  kHz  and  8  kHz.  The 
considered  values  of  the  sea  state  parameters  (a,  T )  are  (0.1333,1.5766)  meters,  which 
corresponds  to  the  wind  speed  of  5  m/sec.  The  considered  dimensions  of  the  patch  are 
X  =  Y  =  10T  ~  15.58  meters. 

The  results  are  plotted  in  the  interval  a\  £  (amin,7r/6)  where  amin  satisfies  the  con¬ 
dition  4 k2a2  sin2  aq  =  gm in  and  gmin  is  set  at  gm ;n  =  2  to  keep  the  magnitude  of  the 
remainder  in  (C14)  small  (see  Figure  Cl  (b)). 

The  plots  of  (101),  (102)  and  (99)  in  Figure  3  (a)  are  reasonably  close  in  the  4  kHz  case, 
whereas  they  become  quite  separated  in  the  8  kHz  case  at  larger  grazing  angles.  Numerical 
integration  of  (75)  is  required  in  order  to  test  whether  or  not  the  dissimilarity  of  plots  in 
Figure  3  (b)  is  caused  by  the  artefacts  of  the  adopted  asymptotic  approximation.  It  also 
would  be  interesting  to  assess  the  difference  between  (67)  and  (74),  and  how  this  difference 
depends  on  the  angle  of  incidence.  Such  tests  are  beyond  the  scope  of  the  current  work 
and  can  be  pursued  in  a  separate  study. 
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Appendix  A  List  of  Notation 


This  section  provides  the  list  of  notation  used  in  Sections  1-5. 


r 

r 

ro 


hi  He 

k 

ksc 

k± 

So 

SM 

Xo,  Vo  and  z0 
® y  and  gz 

no 

-Ro(ro) 
z0  =  h(x0,yo ) 
Oi,  02 ,  03 

A,  B ,  C 
a,  b,  c 

Iff 

sinc(s) 

X,Y 

Am 

kc 

Ic 

Id  and  (Id) 
I  and  (I) 

v ’lc  and  ipsle 

f(o  \,e2,e3) 

p{h) 


the  radius- vector  of  the  observer/receiver  location 
r  =  |  r  | 

the  radius-vector  of  the  surface  integration  variable 

the  complex  amplitude  of  the  incident  pressure  field 

the  complex  amplitude  of  the  scattered  pressure  field 

the  wave  vector  of  the  incident  field 

the  wave  number:  k  =  |kinc|  =  |ksc| 

the  wave  vector  of  the  scattered  field 

these  vectors  are  defined  as  k^  =  kjnc  ±  ksc 

the  scattering  surface,  or  a  patch  of  the  scattering  surface 

projection  of  Sq  on  the  mean  plane 

Cartesian  coordinates  induced  by  the  mean  plane 

unit  vectors  of  the  Cartesian  coordinate  system  induced  by 

the  mean  plane 

a  unit  normal  to  So  at  ro 

the  reflection  coefficient  (here  mostly  assumed  to  be  —1) 
a  Cartesian  representation  of  So 

angles  describing  directions  of  the  incident  and  scattered 
fields,  see  Figure  2 

components  of  vector  k ~ /k,  see  equations  (19)— (21) 
components  of  vector  (k+  —  i?ok_)/l’,  see  equations  (22)— (24) 

the  scattered  field  when  So  is  a  flat  horizontal  rectangular 

plate,  see  equation  (27) 

the  sine  function:  sinc(s)  =  sin(s)/s 

size  parameters  for  rectangular  Sm,  see  (26) 

the  area  of  Sm  in  the  rectangular  case:  Am  =  4 XY 

ipQC  considered  at  specular  scattering  direction  and  Ro  =  —  1, 
see  equation  (29) 

the  intensity  of  the  coherent  component 

the  intensity  and  the  mean  intensity  of  the  diffuse  component 

the  total  intensity  and  its  mean 

the  edge  term  and  the  scattered  field  from  which  the  edge 
term  has  been  removed,  see  (36)  and  (33) 
notation  for  \  (—  +  y  +  c) ,  see  equations  (34),  (37)  and  (54) 
the  probability  density  function  of  the  surface  height  distrib¬ 
ution 
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p2{h1,h2,'R) 

the  probability  density  function  describing  the  two-point 
height  distribution  of  So 

(•) 

the  mean  (expectation)  taken  over  all  permissible  realisations 
of  the  surface  So 

R 

vector  R  =  (Rx,  Ry )  describes  the  displacement  of  the  second 
point  in  the  2-point  height  pdf  with  respect  to  the  first  point 

R 

R=  yjRl  +  R* 

C0(R) 

the  autocorrelation  coefficient,  see  equation  (40) 

T 

the  correlation  length 

X{s)  and  X2(si,s2) 

the  characteristic  functions  associated  with  the  distributions 
p(h)  and  p2(hi,  h2,H),  see  equations  (43)  and  (44) 

X2{si,s2,R) 

same  as  X2(si,  s2) 

Ke( R)  and  P^k) 

autocorrelations  and  power  spectrums  of  the  spatial  processes 
h(x0,yo)  {(■  =  0),  hxo(x0,yo)  =  1)  and  hyo(x0,y0)  (£  =  2) 

h 

I0  =  |^oc(r)|2 

a 

the  standard  deviation  of  surface  heights 

9 

the  roughness  parameter  g  =  k2(T2{cos0\  +  cos 02),  see  equa¬ 
tion  (47) 

90 

the  roughness  parameter  for  the  specular  direction,  that  is 
go  =  4k2a2  cos2  6\ 

Mu) 

the  Bessel  function  of  order  0 

V 

same  as  kr,  v  =  ( vx ,  vy,  vz)  =  k(A,  B ,  C ) 

Vxy 

vxy  =  k\J A2  +  B2,  see  (52) 

S(g) 

notation  for  Jo(vxyT )  e~9^1~expA'r2/T2))  —  e~9  rdr,  an 

integral  factor  in  a  reduced  form  of  {Id),  see  equations  (50) 
and  (61) 

wig) 

the  coefficient  in  the  0{g~ 1)  remainder  of  the  large  roughness 
approximation  of  S(g),  see  equation  (C15) 

Po  and  fi 

angles  associated  with  the  mean  square  wave  slope  and  the 
slope  of  vector  k~,  see  equation  (65)  and  Note  16 

n0 

the  total  rate  of  energy  reflected  from  a  horizontal  flat  rec¬ 
tangular  plate,  see  (67) 

sm(t) 

notation  for  ™  sine2  (mt) ,  see  equation  (68) 

£  and  rj 

auxiliary  variables,  see  equation  (72) 

nd 

diffuse  energy  flux,  see  (75) 

nc 

coherent  energy  flux,  see  (100) 

Q(0  2,  O3) 

auxiliary  function,  see  equations  (77)  and  (78) 

I1 

notation  for  cot2/?o,  see  equation  (79) 

®(02, 03) 

notation  for  (—tan2/?),  see  equation  (80) 

H(02,03) 

the  Hessian  matrix  of  &(02,  O3) 

a\ 

a\  =  tt/2  —  6\ 
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Appendix  B  The  Helmholtz  Integral 


The  Kirchhoff  approximation  is  based  on  representation  of  the  reflected  field  in  terms 
of  the  Helmholtz  integral,  the  derivation  of  which  for  two  different  types  of  the  incident 
field  is  provided  below.  Let  UT( P)  be  the  value  of  a  complex  amplitude  of  a  pressure  field 
at  point  P.  Here  the  subscript  r  is  used  to  indicate  that  UT(P)  will  be  used  to  describe 
the  reflected  component  of  the  total  field.  The  function  UT( P)  satisfies 


V|Ur  +  k2Ur  =  0  ,  (Bl) 

BU 

UT  =  Oir~l)  and  V1  ~  ikUr  =  o(r-1) .  (B2) 

or 

where  r  is  the  distance  from  the  origin.  Equations  (B2)  are  called  the  Sommerfeld  radiation 
condition  [7,  eqn.  (9A-12),  p.  355]. 

The  function 


G( P  -  Q)  =e 


eife|P— Q| 

4vr|P  -  Q 


(B3) 


satisfies  the  equation 


V|G( P  -  Q)  +  k2G{ P  -  Q)  =  <S(P  -  Q) 


(B4) 


and  the  Sommerfeld  radiation  condition. 

Let  S  be  the  boundary  of  a  compact  volume  D  and  consider  C/r(M)  that  satisfies  (Bl) 
outside  D  and  has  the  asymptotic  behaviour  (B2)  as  r  — >  oo.  Let  M  be  a  point  outside 
D,  and  let  T)£)r  be  a  set  comprising  of  a  ball  of  radius  R  centred  at  M,  from  which  we 
have  excluded  the  volume  D  and  an  e-neighbourhood  of  the  point  M.  Also  introduce  the 
notation: 


Br(M),  the  ball  of  radius  r  with  the  centre  at  M, 
S'r(M) ,  the  sphere  of  radius  r  with  the  centre  at  M. 


The  definition  of  V£ir  can  be  written  as  T>£,R  =  Br( M)  \(DU  £L(M)). 

Consider 

by  (B4),  this  is  zero 

r  / — - — ^ — - s 

0  =  lim  /  C/r(P)  [VpG(P  -  M)  +  k2G{ P  -  M)]  dVF  .  (B5) 

Jv,  R 

H — >oo  £'JX 

Use  the  Gauss  (divergence)  theorem  to  transform  the  part  of  the  integral  associated  with 
the  first  term  in  the  square  brackets. 


Z7r(P)VpG(P  -  M)dUP 


VP  (f/r(P)VPG(P  -  M))  -  VPt/r(P)  •  VPG(P  -  M) 


dVp 


UT(P)—G(P-M)dSp  + 

■  OTlp 


d 


Sk(M) 


Ur(P)— G(P-M)dSP 

OTlp 


a 


'  Se  (M) 


Ur(P)~ — G(P  -  M )dSp  -  /  VpUr(P)  •  VPG(P  -  M)dVp 
dn.p  JVeR 


(B6) 


UNCLASSIFIED 


25 


DSTO  -RR-0362 


UNCLASSIFIED 


where  =  np  •  V p  denotes  a  derivative  along  the  outward  normal  np  to  the  integration 
surface  at  point  P. 

Note  that 


lim  <  — 

£ — >0 


Se(M) 


C4(P)w— —  G(P  —  M)  dS-p 

OTlp 


aike  f 


=  lim  /  Ur(P )  — 

47 r 


J5e(M) 

_2^r(M) 


1  i/c 
— o  H - 


d5P 


=  lim  <  47re 


47T 


1  ik 
— £  H - 


=  -C/r(M) 


(B7) 


Next  step 


VPC/r(P)  •  VPG(P  -  MjriUp 


VP  VPUr(P)G(P  -  M)  VpUr(P)G(P  -  M) 


dV p 


(B8) 


Note  1-a.  Since  Vp£/r(P)  +  k2Ur(P)  =  0,  the  second  term  in  the  integrand  of  (B8)  com¬ 
bined  with  the  second  term  in  the  integrand  of  (B5)  gives  zero. 

Use  the  divergence  theorem  to  transform  the  first  term  in  (B8). 

'  9UZ(P) 


>D c 


Vp  (VPUr(P)G(P  -  M))  dV p 
f  dUr(P), 


>SR{  M) 


dn.T 


-G(P-M)  dSP 


'SAM) 


dnp 

dUr(P) 

dnp 


G(P  —  M  )dSP 
G(P  -  M)  dSp  . 


Since  UT( P)  is  regular  at  M  and  |G(P  —  M)|  =  l/(47re)  when  P  G  Pe(M),  we  have 

dUr(  P) 


lim 

£ — >0 


rSf(M) 


dn\ 


-G( P  -  M)  dS-p  =  0  . 


(B9) 


(BIO) 


After  substituting  (B9)  into  (B6),  we  obtain  upon  application  of  the  Sommerfeld  radiation 
condition  that  the  integral  over  S'ft(M)  vanishes  as  R  — >  oo: 


lim 

R—>  oo 


=  lim 

R — »oo 


Sr( m)  L 

r 

Sfl(M) 


[fr(p)Mp-M)_ayp)G(p_M) 


dS* 


dnp  dnp 

—ikUI(P)G(P  -  M)  +  ikUr(P)G(P  -  M)  +  o(R~2)]  dSP 

=  lim  [  [o(R~2)]  dSp  =  0 .  (Bll) 

Hence,  it  follows  from  (B6),  (B7),  (B8),  Note  1-a,  (B9),  (BIO)  and  (Bll)  that  (B5)  is 
equivalent  to 


or 


0  =  -UT(  M)  +  [ 

Js 

Ul{  M)=/r^p 

Js 


-t/,(P)8U^N)  +  8^U0(p-M) 


dnp 

■G(P  -  M)  -  Ur(P 


drip 

dG(  P-M) 
dnp 


dSp 


dSp 


(B12) 
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Note  2-a.  These  derivations  are  essentially  based  on  the  formula 

VpGUr  =  Vp  (vFGUr  -  GVpUr)  +  GVpUr  (B13) 

Substitution  of  (B13)  into  (B5),  taking  account  of  (Bl),  and  application  of  the  divergence 
theorem  to  the  remaining  term  gives 


rj[G,UT]dSp-  /  r][G,Ur}dSP-  /  r][G,  Ur\dSP 
SR{  M)  JSe(  M)  Js 


(B14) 


where 


V[G,  Ur] 


dG{ P  -  M) 
dn-p 


Ur( P)  -  G(P  -  M) 


dUr(P) 

dn-p 


We  have  shown  the  following 


lim  /  r][G,Ur\dSp  =  0  [by  the  Sommerfeld  condition  (B2)] 
R->°°  Jsr(M) 

lim  [  r][G,Ur]dSp  =  Ur(M)  [by  (B7)  and  (BIO)] 

£^°Jse{  M) 

so  (B14)  can  be  written  as 

Ur(M)  =  -  f  n\G,  Ur]dSp 
Js 


(B15) 


(B16) 

(B17) 


(B18) 


which  coincides  with  (B12). 

By  analogy  with  (B5),  consider 

0=  lim  [  G(P  -  Q)  [VpG(P  -  M)  +  k2G(P  -  M)]  dVp,  (B19) 

JgeR 

K— >  oo 

where  Q  /  M,  both  Q  and  M  are  assumed  to  be  outside  D,  but  inside  BR( M),  and 
ge,R  =  Br(  M)  \  (D  U  Se(M)  U  Bs(  Q)). 

Similarly  to  UT( P)  in  T>£,R-,  the  function  G(P  — Q)  satisfies  (Bl)  in  Q£)r  for  any  positive 
e  and  R.  It  also  satisfies  the  Sommerfeld  radiation  condition.  Using  an  argument  similar 
to  that  used  for  (B18)  we  reduce  (B19)  to 


G(M  —  Q)  =  —  /  r]  [G(P  —  M),  G(P  —  Q)]  dSP 

Js 

-lim  [  r)  [G(P  —  M),  G(P  —  Q)]  dSP  ,  (B20) 

£^°Jse(  Q) 

where  r]  is  defined  by  (B15).  The  function  G(P  —  M)  is  regular  at  P  =  Q,  whereas 
G(P  —  Q)  is  singular  (~  |P  —  Q|_1  as  P  — ■>  Q). 


lim  /  r/[G(P-M),G(P-Q)]d5P 

£^°Jse(  Q) 

=  -  lim  [  G{ P  -  M)gG(P  -  Q)  dSp  =  _G(q  _  M) 

£^°Jse{  Q)  dnp 


(B21) 
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Substitute  (B21)  into  (B20)  to  obtain 


(  n  [G(P  -  M),  G(P  -  Q)]  dSp  =  0  . 

Js 


(B22) 


Now  define  U(P)  =  G( P  —  Q)  +  £/r(P)  and  use  equations  (B18)  and  (B22)  to  obtain 


Ur( M)  =  -  [  r,  [G(P  -  M),  U{ P)]  dSp  , 
Js 


(B23) 


or 


Ut(M)=  f 
Js  . 


SU^G(P~M  )-U(F)aG(P-M'> 


dn\ 


drij 


dS] 


p  • 


(B24) 


The  function  G(P  —  Q)  can  be  interpreted  as  an  incident  wave  t/i(P)  from  a  point  source 
at  Q  and,  given  that  the  total  solution  C/(P)  =  £4(P)  +  UT(P)  satisfies  the  boundary 
conditions  at  S,  the  function  £/r(P)  in  this  interpretation  becomes  the  amplitude  of  the 
reflected  wave.  The  right-hand  side  of  equation  (B24)  is  referred  to  (e.g.  [5])  as  the 

Helmholtz  integral. 


Formula  (B24)  also  holds  true  if  we  use  a  plane  wave  for  t/i(P),  that  is 


tfi(P) 


piks  cos  0-\-itp 


(B25) 


where  s  is  the  distance  from  the  M  to  P  and  0  is  the  angle  between  the  wave  vector  k 
and  vector  rjyip  =  tm  —  rp.  For  example,  if  C/i(P)  =  e*k  rp,  then 

—  e*k-(rMp+rM)  _  gifcs cos 0+i(krM) 


The  function  (B25)  doesn’t  have  singular  points  and  satisfies  the  equation  (Bl)  every¬ 
where  in  R3  but  doesn’t  satisfy  the  Sommerfeld  radiation  condition. 

As  in  the  case  P)  =  G(P  —  Q),  we  start  with  an  equation  similar  to  equation  (B19). 


0=  lim  [  C/j(P)  |"VpG(P  —  M)  +  k2G(P  —  M)1  dVp,  (B26) 

JveR 

it— >  OO 

Since  Ui(P)  is  regular  everywhere,  we  use  here  the  same  volume  of  integration  T>£ ^  as  in 
(B5). 

Using  the  same  argument  as  in  the  derivation  of  formula  (B14)  we  obtain 


0=  lim  J  [  ri[G,Ui\dSp-  [  ri[G,Ui}dSF-  [ r)[G,Ui]dSP  \  , 
^_>0  Jsr(  m)  Jsa  m)  Js 


JSe{  M) 

where  we  can  show  as  before  (see  (B17))  that 


(B27) 


lim  /  rj[G,  Ui}dSF  =  UJM)  =  . 

e^°  Js£{ M) 


(B28) 


Because  U;(P)  doesn’t  satisfy  the  Sommerfeld  radiation  condition  the  integral  over  S'r(M) 
doesn’t  vanish  as  R  —>  oo.  Let  us  show  that 


lim  f  rj[G,Ui\dSP  =  eiip .  (B29) 

R->0°  JsR( M) 
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Since,  on  S'r(M), 


/  piks  \ 

iks  cos  Q-\-itp  _  I  _ _ _  |  w  Q^ks  cos  9-\-i<p 


g iks  cos  6-\-iip-\-iks 

An 


0iks  ' 

4-rrs  )  ds ' 


—  -  4  -  -ifccosfl 
s  sz  s 


s—R 


s—R 


we  have 


isR(  m) 


r}[G,UddSP  = - 


2TrR2eikR+iv  fn 
47 tR  Jn 


d9  sin  9eikRcos  8 


ik  —  —  —  ik  cos  9 
R 


1  r1 

=  -RelkR+icP  /  d£eikRt 

2  J- i 


ik£  +  —  —  ik 

K 


ikR£  pikR £  pikR £  pikR£ 


R 


ikR 2  ifci?2  /? 


«=i 


«-•  i 


_  _  r  ikR  —ikR  _  ikR  .  —  _  icp 

D  C  1  C  D  I  D  C  < 


which  proves  (B29). 

Substitution  of  (B28)  and  (B29)  into  (B27)  gives 


/  rj[G,  Ui]dSp  =  0 .  (B30) 

Jse(  M) 

Combining  (B30)  with  (B18)  and  using  U(P)  =  £/i(P)  +  Ur( P)  we  obtain  that  (B23)  and 
(B24)  are  also  valid  for  E/i(P)  =  eifa,cos®+i¥\ 
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Appendix  C  Behaviour  of  S(g)  at  large  g 


In  this  section  we  estimate  the  behaviour  of  (61)  at  large  values  of  the  roughness 
parameter  g.  We  obtain  both  the  principal  term  of  the  asymptotic  expansion  and  an 
estimate  of  the  remainder. 

Without  loss  of  generality  we  can  assume  that  T 2  =  Pg ,  where  P  =  P(g)  is  a  positive 
parameter.  Equation  (61)  is  equivalent  to 


S(g) 


' e~r2/P  +  e-g{l-eM-P/T^))  _e-P/P  _e-g 


rdr . 


(Cl) 


Because 

Jo{vxyT)e~T2/PTdT  =  ^Pexp  (-^v2xyP^j 
we  obtain  that  (Cl)  takes  the  form 


S(g) 


-P  exp 


v2  P 

uxy± 


+  H(g ) 


(C2) 


where 


H(g)  =  /  Jo{vxyT ) 


e-g(l-exp(-r2/T2))  _  &-r 2/P  _  &-g 


rdr . 


We  will  show  that,  when  g  — ■>  oo, 


\H{g)\  <  —P  (^  H(g)  =  O^P) . 

g 


Use  |  Jo(i)  |  <  1  to  write 

l^(5)|  < 


,-ff(l-exp(-r2/Pg))  _  -t2/P  _  -g 


rdr . 


(C3) 


(C4) 


In  (C4)  change  the  integration  variable  from  r  to  £  =  ge  T~/Pg ,  for  which  we  have 

=>  rdr  =  --di 

and 

£(0)=fo  £(°o)  =  0  . 

As  a  result,  the  relation  (C4)  is  transformed  to 

'O9 


“ft 


,-<?+£ 


-  I  -  I  -  e 

g 


-g 


dt, 


AS'-- 


(C5) 


30 


UNCLASSIFIED 


UNCLASSIFIED 


DSTO-RR-0362 


where  we  assumed  that  g  >  1  and  have  used  the  fact  that,  for  £  £  [0 ,  <7] , 


>  g— 


4A  e' 


<>i*Y  «  « 


-9+C  >  f  £ 
9 


Z'9 


9  V  9 

Rewrite  (C5)  as 

I#  (<?)|  <  #1(5)  +  ^2(5)  +  +  ^4(5) 

and  estimate  each  term  on  the  right  separately. 

Hi(»)  =  <  ^~>(e  -  1) 

=  H  ^e~9H<*  "  f  (X  -  S~fl)  =  ^3(5)  +  ^3b(ff)  • 

Use  integration  by  parts  to  estimate  H^{g): 

Pge~ 


Hl(9)  = 


-9 

e € 

9  e*| 

9  r eC 

~ 

I 

+  72 

1  S 

pi  ¥(. 

P 

2 


1  —  ge  9+1  +  -  —  <?e  9+1  +  25c  9 
9 


ra  ei 


K  ed\ 


Next  estimate  the  last  term  in  (CIO): 

r9  Pi  r9/2  gff/2 


ge 


-9 


h  e 


-pid£  <ge  9 


e 


d£  +  ge 


-9 


r9  ei 


4/2  (9/ 2): 


=  9C  9/2  (  \  ~  -g2  )  +  ^  f1  “  6  9/2)>  =  ^ 


Substitution  of  (Cll)  into  (CIO)  gives 

H3(g)  =  ^[l  +  0(g-1)}- 
Estimate  the  last  term  in  (C6): 

r9  Pge-9 


Hi{g)  = 


1 1 


2£ 


Combining  (C6),  (C7),  (C8),  (C12)  and  (C13),  and  using  the  fact  that 

ge~9,  g~9,  glog(g)e~9 

vanish  faster  than  0(g~n),  where  n  is  an  arbitrarily  chosen  number,  we 
Substitute  (C3)  into  (C2)  and  use  the  fact  that  P  =  T2 / g  to  obtain 


(C6) 

(C7) 

(C8) 

(C9) 


(CIO) 


(Cll) 


[l-g~9]  =  0(g~1)P . 

(C12) 

Pge~ 9  , 

=  2  log  5- 

(C13) 

T2 

Sis)  ~  2? 


exp 


y2  T 2 

xy 


4  g 


+  0(g 


-u 


obtain  (C3). 

(C14) 


as  g  —y  00. 
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Figure  Cl:  Comparison  of  the  exact  and  asymptotic  forms  of  S(g)  for  the  specular  case 
and  plots  of  the  remainder  w(g)/g  (b)  and  its  coefficient  w(g)  (c). 


Note  20.  The  remainder  0{g  *)  in  (C14)  is 

1 0{g~l)\  =  w(g)/g , 


(C15) 


where  w(g)  is  a  bounded  function.  We  have  verified,  using  (C6)  and  numeric  evaluation 
of  (C7)-(C11)  and  (C13),  that  \w(g)\  <  6.2  for  all  g  >  1,  |u;(g)|  becoming  less  than  1.1 
when  g  >  200.  These  are  quite  coarse  estimates,  in  reality  w(g)  can  be  smaller.  A  more 
accurate  estimate  for  w(g)  is  obtained  by  considering 


w{g)  =  g 


2 g_ 

rp2 


S(g)  -  exp 


vlyT'2\ 

) 


where  S(g)  is  evaluated  numerically  using  either  (61)  or 


S(g) 


m=l 


gm  exp 


v2  t2\ 

uxy~L  \ 

4m  I 


(see  [6,  eqn.  (35),  p.  86]).  Figure  Cl  compares  the  exact  and  asymptotic  forms  of  S(g)  for 
the  specular  case.  Figure  Cl  (a)  shows  that  S(g)  and  its  asymptotic  approximation  are 
quite  close  even  at  moderately  small  (~  1)  values  of  g.  The  plot  of  the  remainder  (C15) 
in  Figure  Cl  (b)  shows  that  it  does  not  exceed  0.5  at  g  >  1.  Figure  Cl  (c)  confirms  our 
upper  estimates  for  w(g)  and  shows  that  for  the  specular  case  w(g)  <  1.5  if  g  >  1  and 
w(g)  <  1.05  if  g>  50. 
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